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\^ ■ The physics of compact objects is one of the very actively devel- 

, oping branches of theoretical investigations, since the careful analysis 

of different models for the internal structure of these objects, their 
evolutionary characteristics and the comparison with modern obser- 
Q^i vational data give access to discriminate among various speculations 

^ I about the state of matter under extreme conditions. 

' The lecture provides an overview to the problem of neutron star 

c/3 . cooling evolution. We discuss the scheme and necessary inputs for neu- 

^ ' tron star cooling simulations from the background of a microscopic 

modeling in order to present the surface temperature - age charac- 
^ ■ teristics of hybrid neutron stars including the determination of main 

^ ! regulators of the cooling process, the question of neutrino production 

and diffusion. 
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1 Introduction 



The interiors of compact stars are considered as systems where high-density 
phases of strongly interacting matter do occur in nature, see Shapiro and 
Teukolsky P], Glendenning |2] and Weber |3j for textbooks. The conse- 
quences of different phase transition scenarios for the coohng behaviour of 
compact stars have been reviewed recently in comparison with existing X-ray 
data |4j. 

The Einstein Observatory was the first that started the experimental 
study of surface temperatures of isolated neutron stars (NS) . Upper limits for 
some sources have been found. Then ROSAT offered first detections of sur- 
face temperatures. Next X-ray data came from Chandra and XMM/Newton. 
Appropriate references to the modern data can be found in recent works by 
[HI ini HI E] ; devoted to the analysis of the new data. More upper limits and 
detections are expected from satellites planned to be sent in the nearest fu- 
ture. In general, the data can be separated in three groups. Some data show 
very "slow cooling" of objects, other demonstrate an "intermediate cooling" 
and some show very "rapid cooling". Now we are at the position to carefully 
compare the data with existing cooling calculations. 

The "standard" scenario of neutron star cooling is based on the main 
process responsible for the cooling, which is the modified Urea process (MU) 
nn —>■ npeu calculated using the free one pion exchange between nucleons, see 
Pj. However, this scenario explains only the group of slow cooling data. To 
explain a group of rapid cooling data "standard" scenario was supplemented 
by one of the so called "exotic" processes either with pion condensate, or 
with kaon condensate, or with hyperons, or involving the direct Urea (DU) 
reactions, see fOl |Tj and refs therein. All these processes may occur only 
for the density higher than a critical density, (2 6) uq, depending on the 
model, where no is the nuclear saturation density. An other alternative to 
"exotic" processes is the DU process on quarks related to the phase transition 
to quark matter. 

Particularly the studies of cooling evolution of compact objects can give 
an opportunity for understanding of properties of cold quark gluon plasma. 
In dense quark matter at temperatures below ~ 50 MeV, due to attractive in- 
teraction channels, the Cooper pairing instability is expected to occur which 
should lead to a variety of possible quark pair condensates corresponding to 
color superconductivity (CSC) phases, see [TJj for a review. 

Since it is difficult to provide low enough temperatures for CSC phases in 
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heavy-ion collisions, only precursor phenomena ^21 ^] are expected under 
these conditions. 

CSC phases may occur in neutron star interiors ^3] and could manifest 
themselves, e.g., in the cooling behavior fTEl IT^ IT71 ITH]. 

However, the domain of the QCD phase diagram where neutron star 
conditions are met is not yet accessible to Lattice QCD studies and theoret- 
ical approaches have to rely on non-perturbative QCD modeling. The class 
of models closest to QCD are Dyson-Schwinger equation (DSE) approaches 
which have been extended recently to finite temperatures and densities P^ . 
Within simple, infrared- dominant DSE models early studies of quark stars 
pUj and diquark condensation [21j have been performed. 

Estimates of the cooling evolution have been performed jT^] for a self- 
bound isothermal quark core neutron star (QCNS) which has a crust but no 
hadron shell, and for a quark star (QS) which has neither crust nor hadron 
shell. It has been shown there in the case of the 2SC (3SC) phase of QCNS 
that the consequences of the occurrence of gaps for the cooling curves are 
similar to the case of usual hadronic neutron stars (enhanced cooling). How- 
ever, for the CFL case it has been shown that the cooling is extremely fast 
since the drop in the specific heat of superconducting quark matter dom- 
inates over the reduction of the neutrino emissivity. As has been pointed 
out there, the abnormal rate of the temperature drop is the consequence of 
the approximation of homogeneous temperature profiles the applicability of 
which should be limited by the heat transport effects. Page et al. (2000) ^Hl 
estimated the cooling of hybrid neutron stars (HNS) where heat transport ef- 
fects within the superconducting quark core have been disregarded. Neutrino 
mean free path in color superconducting quark matter have been discussed 
in 1221 where a short period of cooling delay at the onset of color supercon- 
ductivity for a QS has been conjectured in accordance with the estimates of 
|,15|j in the CFL case for small gaps. 

A completely new situation might arise if the scenarios suggested for 
(color) superconductivity [23 121| besides of bigger pairing gaps (A^ ~ 50 
100 MeV) will allow also small diquark pairing gaps (A^ < 1 MeV) in quark 
matter. 

The questions which should be considered within these models are the 
following: (i) Is strange quark matter relevant for structure and evolution 
of compact stars? (ii) Are stable hybrid stars with quark matter interior 
possible? (iii) What can we learn about possible CSC phases from neutron 
star cooling? Further on in this lectures we discuss the scheme and the results 
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of realization of the these points in relation with the cooling evolution of 
compact objects. 

In the consideration of the scenario for the thermal evolution of NS and 
HNS we include the heat transport in both the quark and the hadronic 
matter. We will demonstrate the influence of the diquark pairing gaps and 
the hadronic gaps on the evolution of the surface temperature. 

The main strategy of the simulation of the cooling evolution of compact 
objects is presented in Fig ^ On the top of scheme we have the general 
theoretical background of QCD models as it has been discussed in the intro- 
duction. On the second level of the scheme we separate two branches one for 
the structure of the compact objects and the other for the thermal proper- 
ties of the stellar matter. On the bottom of the scheme those two branches 
are combined in the code of the cooling simulations project and entail the 
comparison of the theoretical and observational results. 



Microscopic Approach to 
Dense QCD / Nuclear Matter 



EoS and NS constraints 



NS configurations 



Pairing Gaps 



Neutrino Emissivities, ^ 
Transport Properties, 



Cooling Code 



Obs. Data 



Figure 1: The scheme linked problems and tasks, which has to be done to 
simulate the cooling of compact objects 
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2 EoS of stellar matter and Compact Star 
configurations 

2.1 Hadronic matter EoS and DU threshold 

For the modeling of the dense hadronic matter different approaches, hke 
relativistic and non-relativistic Briickner-Hartree-Fock approach including 
the three body correlations |25^ EEl I2Z1 12H] and different parameterizations 
of relativistic mean field models, are discussed successfully in the literature 

For the appropriate cooling simulations we exploit the EoS of (specif- 
ically the Argonne V18 + 6v + UIX* model), which is based on the most 
recent models for the nucleon-nucleon interaction including a parameterized 
three-body force and relativistic boost corrections. Actually we adopt a sim- 
ple analytic parameterization of this model by Heiselberg and Hjorth- Jensen 
[2H1, hereafter HHJ. 

The latter uses the compressional part with the compressibility K ~ 
240 MeV, a symmetry energy fitted to the data around nuclear saturation 
density and smoothly incorporates causality at high densities. The density 
dependence of the symmetry energy is very important since it determines 
the value of the threshold density for the DU process. The HHJ EoS fits the 
symmetry energy to the original Argonne V18 + 6v + UIX* model yielding 
~ 5.0 no (Mf^ ~ 1.839 M©). 

Fig. 121 demonstrates the partial densities of p, e and fi^ in HHJ and 
relativistic non-linear Walecka (NLW) model in the parameterization of |H2] , 
adjusted to the following bulk parameters of the nuclear matter at saturation: 
no = 0.16 fm~^. Anyhow, in the given NLW model the threshold density for 
the DU process is n^^ ^ 2.7 no . 

2.2 Quark matter EoS with 2SC superconductivity 

The quark matter models EHl iSlj predict, that the diquark pairing con- 
densate is possible for temperatures and densities relevant for compact ob- 
jects. The order of magnitude of pairing gaps is 100 MeV and a remarkably 
rich phase structure of the matter has been identified [TTl ESI EZI- 

This leads to expectations of observable effects of color superconducting 
phases, particularly in the compact star cooling behaviour [13 CHI 113 UH 
EHl- Generally color superconductivity is involved in all aspects of neutron 
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Figure 2: Densities of the charged particles as a function of baryon density 
for HHJ and NLW model EoS. Possibility of pion condensation is suppressed. 



star studies, such as magnetic field evolution HOI HH |12] or burst-type 
phenomena |13 EH 113 EH] • 

The applications of the nonlocal chiral quark model developed in jl7j for 
the case of neutron star constraints shows that the relevant CSC phase is a 
2SC phase while the omission of the strange quark flavor is justified by the 
fact that chemical potentials in central parts of the stars do barely reach the 
threshold value at which the mass gap for strange quarks breaks down and 
they appear in the system jJSl US] • 

It has been shown in that work that the Gaussian formfactor ansatz of 
quark interaction (hereafter we call it SM model) leads to an early onset of 
the deconfinement transition and such a model is therefore suitable to discuss 
hybrid stars with large quark matter cores [SU] . 

The resulting quark matter EoS within this nonlocal chiral model |17] 
can be represented in a form reminiscent of a bag model 

P^'^ = Mi^b) - B('\fiB) , (1) 

where PidifJ^B) is the ideal gas pressure of quarks and B^^\iib) a density 
dependent bag pressure, see Fig. El The occurrence of diquark condensation 
depends on the value of the ratio rj = G2/G1 of coupling constants and 
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the superscript s E {S, N} indicating whether we consider the matter in 
the superconducting mixed phase (r^ = 1) or in the normal phase [r] = 0), 
respectively. 




Figure 3: Bag pressure for different formfactors of the quark interaction in 
dependence on the baryon chemical potential for rj = and for rj = 1. For 
the latter the superconducting phase is realized. 



2.3 EoS of hybrid star matter 

In some density interval above the onset of the first order phase transition, 
there may appear a mixed phase region between the hadronic (confined quark 
phase) and deconfined quark phases, see ^1]. Ref. disregarded finite 
size effects, such as surface tension and charge screening. On the exam- 
ple of the hadron-quark mixed phase Refs jH2] demonstrated that finite size 
effects might play a crucial role substantially narrowing the region of the 
mixed phase or even forbidding its appearance in a hybrid star configura- 
tion. Therefore we omit the possibility of the hadron-quark mixed phase in 
our model assuming that the quark phase arises by the Maxwell construction. 

To demonstrate the EoS with quark-hadron phase transition applicable 
for hybrid star configurations in the next Section, we show in Fig. H] re- 
sults using the relativistic mean field (RMF) model of asymmetric nuclear 
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matter including a non-linear scalar field potential and the p meson (non- 
linear Walecka model) ( see j2]) for the cases 77 = 1 (left panel) when the 
quark matter phase is superconducting and for 77 = (right panel) when it 
is normal. 




500 1000 1500 500 1000 1500 
E [MeV/fm^] e [MeV/fm^] 

Figure 4: EoS for strongly interacting matter at zero temperature under 
compact star constraints for the coupling parameter 77= 1 (left panel) and 
?7 = (right panel). Dashed line: relativistic mean- field model for hadronic 
matter; solid, dash-dotted and dotted lines correspond to quark matter with 
Gaussian, Lorentzian and NJL formfactor functions, respectively. 



In case the hadonic EoS is chosen to be HHJ model and the quark one 
is SM model we found a tiny density jump at the phase-boundary from 
^hadr ^ fm~'^ to n'^"'^^^ ~ 0.46 fm~^. The critical mass, when the branch 
of stable hybrid stars starts is M^^"^'^'^ = 1.214 M©. 

As shown in Refs. |53| 1321 luark matter itself can exist in a mixed phase 
of 2SC and normal states. The presence or absence of the 2SC - normal 
quark mixed phase instead of only one of those phases is not so important 
for the hybrid star cooling problem since the latter is governed by processes 
involving either normal excitations or excitations with the smallest gap. 
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2.4 Stability of Hybrid star configurations 



Not all compositions of hadronic and quark matter EoS's give stable configu- 
rations for compact objects, which could be demonstrated using the HHJ 
model EoS for hadronic and the two flavor nonlocal chiral model EoS for 
quark matter. 

The mechanical structure of the spherically symmetric, static gravita- 
tional self-bound configurations of dense matter can be calculated with the 
well-known Tolman-Oppenheimer-Volkoff equations. These are the condi- 
tions of hydrostatic equilibrium of self- gravitating matter, see also 



Here e{r) is the energy density and P{r) the pressure at distance r from 
the center of the star. The mass enclosed in a sphere with radius r is defined 



These equations are solved for given central baryon number densities, 
nsir = 0), thereby defining a sequence of star configurations. 

The configurations are stable if they are on the rising branch of the mass 
- central density relation. 

The relation between pressure and energy density is given by the choice 
of the corresponding EoS model, which generally can be a function of tem- 
perature. The temperature is either constant for isothermal configurations, 
see PP, or has some profile, when the cooling evolution is assumed to be a 
hydrodynamically quasi-stationary process. 

Hot quark stars have been discussed, e.g., in fIU\ ITTl 13^ . 

For late cooling, when the temperatures are below T < 1 MeV the effects 
of temperature on the distribution of matter are negligible and the calcula- 
tions for the star structure can be performed for the T = case. 

Dots in Fig. |21 indicate threshold densities for the DU process. The 
possibility of charged pion condensation is suppressed. Otherwise for n > 
nj?^ the isotopic composition is changed in favor of increasing proton fraction 
and a smaller critical density for the DU reaction. Deviations in the M{n) 
relation for HHJ and NLW EoS are minor, whereas the DU thresholds are 
quite distinct. 



dP{r) 



dr 



[e{r) + P{r)] [m{r) + 4TTr^P{r)] 
r[r — 2m(r)] 



(2) 



by 




(3) 
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In Fig. we show the mass-radius relation for hybrid stars with HHJ 
EoS vs. Gaussian nonlocal chiral quark separable model (SM) EoS. Config- 
urations with possible 2SC phase given by the solid line, are stable, whereas 
the hybrid star configurations without 2SC large gap color superconductivity 
(dash-dotted lines) are unstable. In case of "HHJ-SM with 2SC" the max- 
imum neutron star mass proves to be 1.793 Mq. For an illustration of the 
constraints on the mass-radius relation which can be derived from compact 
star observations we show in Fig. El the compactness limits from the ther- 
mal emission of the isolated neutron star RX J1856. 5-3754 as given in 
and from the redshifted absorption lines in the X-ray burst spectra of EXO 
0748-676 given in j^. These are, however, rather weak constraints. 

2 
1.8 
1.6 
^ 1-4 

1 

0.8 
0.6 

9 10 11 12 13 

R [km] 

Figure 5: Mass - radius relations for compact star configurations with dif- 
ferent EoS: purely hadronic star with HHJ EoS (dashed line), stable hybrid 
stars with HHJ - Gaussian nonlocal chiral quark separable model (SM) with 
2SC phase (solid fine) and with HHJ - SM, without 2SC phase (dash-dotted 
line). Data for two sources are also indicated, see IHBj. 
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3 Heat transport, Neutrino production and 
diffusion 



3.1 Boltzmann equation in curved space - time 

For study of cooling evolution, we will use the diffusion approximation for 
neutrino transport, because it allows us easily to assess how global character- 
istics of the neutrino emission such as average energies, integrated fluxes, and 
time scales change in response to various input parameters. This approach is 
sufficient to establish connections between the microphysical ingredients and 
the duration of deleptonization and cooling time scales, which are necessary 
to estimate possible effects on neutrino signals. However, the detailed results 
of simulations for early cooling will not be presented in this lecture. The 
reader is referred to the works jHZi EH] and citations within. 

We will assume that the neutron star has an interior where the matter is 
spherically symmetric distributed. The evolution of the star will be assumed 
to be quasi stationary, which means that for each time step the velocity of 
the matter is zero. 

The Boltzmann equation (BE) for massless particles is 



coll 



where / is the invariant neutrino distribution function, p"" is the neutrino 
4-momentum and F^^ are the Christoffel symbols for the metric 

dr^ = -e^Ue + e^^dr^ + r^dO^ + sin^ 6 (5) 

of static spherically symmetric space-time manifold. 

For simplicity one can use the comoving basis and rewrite the BE 
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where are basis vectors in a comoving frame to an observer and in the 
static case they are diagonal. Indices a, 6, c are running from to 3. The F^^ 
are Ricci rotational coefficients with the non-zero components 

T^l _ -PO _„-A^ /yN 

J- 00 — 01 — ) \' ) 
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p2 _ _pl _ _pl _ p3 _ 

-■-21 ~ 22 ~ 33 ~ -"^ 31 ~ ^ ■ 

33 ~ 32 ~ 



r 

The neutrino 4-momentum is 

= (u, LUjj,, uj{l - /i^)^^^ cos uj{l - /i^)^^^ sin $ 

where n is the cosine of the angle between the neutrino momentum and the 
radial direction, u is the neutrino energy in a comoving frame. Using the 
previous definitions we will have 

-eo^ + -/^^i9;:--/^roo^--(i-^)(roo + r22)^=(^^jj8) 



3.1.1 Equations in Linear response Approximation 

After applying the operator 

1 

-J^ dfifi\ « = 0, 1,2, 
to equation (jHJ and defining the i^^ moments 
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we get for i = 



^ (^0^ + ^^^) - I ^00^ 1 - 2-(r^o + rL)Mi = Qo, (9) 

and for i = 1 



1^ <^Q—. — \- e 



- (^00^) + ^(r^o + rL)(Mo - 3M2) = Q, .(10) 



Let us introduce Ny, F^,, and Sat as the number density, number flux and 
the number source term, respectively. 



i ^AW, K = l ^M.^^ s„ = /^ ^o„.. (11) 
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while J^^ Hy, Py, and 5"^; are the neutrino energy density, energy flux, pres- 
sure, and the energy source term 

JO 27r^ Jo /vr^ 

P. = r^,M,u\ S,= r^^Qou;\ (12) 

Jo ZTT^ JO ZTT^ 

After integration over the neutrino energy and utilizing the continuity equa- 
tion under assumption of a quasi-static evolution, one can recover the well 
known neutrino transport equation ^60j 



dt da hb' 

dt " dt da He 

The distribution function in the diffusion approximation could be repre- 
sented in the following way 

/(a;,/x) = /oM+/i/iM, /o = [l + e(^)]-S (15) 

where /o(i-i-') is the distribution function in equilibrium (T = T^at, IJ^v = 
/i^*^) with the neutrino energy uj and the chemical potential respectively. 
Hereafter the dependence of /o and /i on oj and non explicit dependence on 
space-time coordinates will be assumed without listing in arguments. 
Thus the moments Mj of the distribution function / are 

Mo = /o, Mi = ^/i, M2 = ^/o, and M3 = ^/i . (16) 
Therefore the equation (fTUj) now reads 

\dr dr duj ) uj 
The collision term Qi can be represented as 
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Here ja is the emissivity, is the absorptivity, js and are the scattering 
contributions. 
Namely 

1 /"OO rl r2n 

(27r) "'0 J-i Jo 

1 1 rOD rl /•2lT 

T = 7:^ ^^'^'} ^-"7 d^[l-f{u\^i')]RT\u:,uJ\cose), (20) 

As (27r) "'0 J-i Jo 

where 6 is the scattering angle. The relation between emissivities and ab- 
sorptivities (i?*" and are scattering kernels) is given by 

= e^("-'^^')j„(cu) and = e^(^'-^)i?r* • (21) 

Aa(t^) 

Using the Legendre expansion for the moments one has 

= J^^ d cos e P;(C0S e)R'^\uJ, U', COS 9) . (22) 

Performing the angular integrations of Eq. ()18p one can define the relation 
between Qo, Qi and i?f *. 

After substitution of the expression for Qi into Eq. lfTTj) the relation be- 
tween /o and /i reads 



dfo _ _d(j) dfo 
dr dr duj 



D{uj) 



Ja + ^ + t^l 



(23) 



where D{ijj) is the diffusion coefficient. Using the relation between ^ and 
^ and the notation 77 = /ii./T to be the neutrino degeneracy parameter, we 
obtain 



-D(uj)e 



T 



dr] uj d{Te^ 



dr Te't' dr 
Now the energy-integrated lepton and energy fluxes are 



dfo 
duj 



(24) 



67r2 



67r2 



d{Te^) 

dr 
d{Te^) 



drj 
dr 



dr ~^ ^ ^ ^ dr 



(25) 
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Here the coefficients D2, D3, and D4 are defined through the integrals over 
^as , , /-oo 







The exphcit expression for the absorption mean free path is 

Xa T^'^ Jo 



/ dE, El [1 - / cos ^ + i?2 + 5i?3](26) 

Jo J-1 i — e ^ 



i?-* = (cos g 1) ^ ^ 
1 — e ^ 

The response functions written in terms of the polarization functions are 

= (V2 + ^2) [j^nf(go,g) + Imn^(go,g)] (28) 
i?2 = (V' + ^2) Imn^(go,g)-^' Imn^(go,g) (29) 
R, = 2V^Imn^^(go,g), (30) 

where the coupling constants for the absorption are V = gy cos(6'c) and A = 
qa cos(6'c) and for the scattering V = (y'y/2 and A = Qa/'^- is the Cabibbo 
angle, gy and qa are the vector and axial-vector couplings. 

Introducing the s entropy per baryon and Yl = {Ni, + Ne + N^)/nB lepton 
number fraction and with help of conservation laws for energy and baryon 
number combined with Eqs. (fT^ and (fT^ we finally will obtain the energy 
transport equation 



and the equation for lepton diffusion 

dt da 

The energy and lepton F^ fluxes are defined in Eq. ip^j) . 

These coupled equations for the unknowns T and n^, with the Eq.Q 
for the star structure completely describe the early cooling evolution of the 
compact objects in quasi-stationary regime. The solutions and more detailed 
information one can find in [2Z]. In the next section we will use the energy 
transport equation for consideration of late cooling evolution when the neu- 
trinos are already untrapped, i.e. the mean free path is larger than the 
star radius R and fi^, = 0. 
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3.2 Late cooling evolution 

The cooling scenario is described with a set of coohng regulators. In the 
regime without neutrino diffusion in the righthand side of Eq. (fT^ and 
therefore also in Eq. for energy transport we introduce the energy loss 
term due to neutrino emission additional to Se and instead of the Hy we 
use the total energy flux function l{r,t). The temperature profile T{r,t) is 
the one time dependent unknown dynamical quantity and during the cooling 
process due to an inhomogeneous distribution of the matter inside the star 
and the finite heat conductivity it can differ from the isothermal one. The 
heat conductivity k (the corresponding term to D4 in Eq. (jH^) determines 
the relation between l{r,t) and gradient of T{r,t). In isothermal age the 
temperature on the inner crust boundary (T^) and the central temperature 
(Tc) are connected by the relation Tm = Tcexp[(f){0)—(f){R)]. Here 0(0)— 0(i?) 
is the difference of the gravitational potentials in the center and at the surface 
of the star, respectively. 

The flux of energy l{r) per unit time through a spherical slice at the 
distance r from the center, is proportional to the gradient of the temperature 
on both sides of this slice, 



l{r) = -47rr'K{r)^-^e-'^\l 1 - — , (33) 



where the factor e'^^y 1 — ^ corresponds to the relativistic correction of the 
time scale and the unit of thickness. The equations for energy balance and 
thermal energy transport are [3] 

A(/e-) ^ -l(e.e- + c4(Te^)), (34) 

where n = n{r) is the baryon number density, A = A{r) is the total baryon 
number within a sphere of radius r. One has 



dr 1 2M 



OA Anr'^n V r 

The total neutrino emissivity e^, and the total specific heat cy are given as the 
sum of the corresponding partial contributions defined in the next subsec- 
tions. The density profiles ni{r) of the constituents i of the matter are under 
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the conditions of the actual temperature profile T{r,t). The accumulated 
mass M = M{r) and the gravitational potential = 0(r) can be determined 
by 



dM 
'dl 




Airry + M 



OA A-nr'^n . Ii_m. 



(37) 
(38) 



where the energy density profile e = e{r) and the pressure profile p = p{r) is 
defined by the condition of hydrodynamical equilibrium (see Eq. ^ 

The boundary conditions for the solution of ()34|) and (jH^j) read l{r = 0) = 
l(A = 0) = and T{A{r^) = A,t) = T{r^ = R^t) = T^{t), respectively. 

In our examples we choose the initial temperature to be 1 MeV. This 
is a typical value for the temperature Topadty at which the star becomes 
transparent for neutrinos. Simplifying we disregard the neutrino influence 
on transport. These effects dominate for t < 1 100 min, when the star 
cools down to T < Topadty and become unimportant for later times. 



3.2.1 Neutrino Processes in dense matter 

We compute the NS thermal evolution adopting our fully general relativistic 
evolutionary code. This code was originally constructed for the description 
of hybrid stars by ^\ . The main cooling regulators are the thermal conduc- 
tivity, the heat capacity and the emissivity. 



Emissivity 

The luminosities are calculated using the corresponding matrix element of 
the neutrino production process 

(40) 

The most effective process of neutrino production is the direct Urea (DU) 
one, n — > p + e~ + Ve-, which is allowed when the momentum conservation 
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holds pF^n = + PF,e "^^ |PF,n| < |PF,p| + |PF,e|- Siiice DU piocess is 
so intensive that as soon as it takes place the star cools too fast for better 
correspondence between the cooling simulations and the observational data, 
it is likely to assumed DU to be forbidden in the interior of the star. 

The other processes, i.e. the modified Urea (MU), pair breaking and 
formation (PBF) and bremsstrahlung are the main ones governing the cooling 
in hadronic matter. The modified Urea (MU) nn —>■ npeu, np —>■ ppeu, has to 
be medium modified Urea (MMU) process due to the softening of in-medium 
pion propagator. The final emissivity is given by 122] 

e-^^u ^ 8.6 X 102im:tiu(>^ew)'^'CnMuT9Wm-'s-\ (41) 
eP^u = 8.5 X 10'^m;t,^iY,uy%MvT^ergcm-h-\ (42) 

Here m* = \Jm*l-^ ^ + pp ^ is the non-relativistic quasiparticle effective mass 
related to the in-medium one-particle energies from a given relativistic mean 
field model for i = n,p. We have introduced the abbreviations rn'^jj = 
{m^/mnY^mp/mp) and m*^^^ = {mymp)^{m*^/mn)- The suppression fac- 
tors are CnMU = CnCp - exp{-[A„(T) + Ap(T)]/T}, (pMU - Cp, and should be 
replaced by unity for T > Tcrn^i, when for given species i the corresponding 
gap vanishes. For neutron and proton S-wave pairing is Aj(0) = 1.76 Tcnt,i 
and for the P-wave pairing of neutrons A„(0) = 1.19 Tcrit,n (see Fig. Ej). 
The function of temperature is given by the interpolation formula 

A^(T) = A(0)y^l - r/T„it,7v. 

To be conservative we have used in ()4H1 the free one-pion exchange esti- 
mate of the NN interaction amplitude. Restricting ourselves to a qualitative 
analysis we use here simplified exponential suppression factors Q. In a more 
detailed analysis these ^j-factors have prefactors with rather strong temper- 
ature dependences [221 • At temperatures T ~ their inclusion only slightly 
affects the resulting cooling curves. For T ^ the MU process gives in 
any case a negligible contribution to the total emissivity and thereby corre- 
sponding modifications can again be omitted. Also for the sake of simplicity 
the general possibility of a ^P2(|'t?.j| = 2) pairing which may result in a 
power-law behaviour of the specific heat and the emissivity of the MU pro- 
cess inn EH inS] is disregarded since mechanisms of this type of pairing are up 
to now not elaborated. Even more essential modifications of the MU rates 
may presumably come from in-medium effects which could result in extra 
prefactors of 10^ -h 10'^ already at T ~ Tg. 
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Figure 6: Neutron and proton pairing gaps according to model I (thick solid, 
dashed and dotted lines) and according to model II (thin lines), see text. The 
ISo neutron gap is the same in both models, taken from 



In order to estimate the role of the in-medium effects in the A^A^ interac- 
tion for the HNS cooling we have also performed calculations for the so-called 
medium modified Urea (MMU) process [HSl IMI by multiplying the rates 1)411] 
by the appropriate prefactor 

^MMU/^MU ^ ^q3 [r6(^')/^8(^ ^ ^^^j ^10/3^ (43) 

where the value T{g') ~ 1/[1-|-1.4 u^^^] is due to the dressing of itNN vertices 
and uj < is the effective pion gap which we took as function of density 
from Fig. 2 of [HZj- 

For T < Tcrit the most important contribution comes from the neutron 
[nSl in2] and the proton [^21 pair breaking and formation processes. We take 
their emissivities from Ref. [UZj, which is applicable for both cases, S- and 
P-wave nucleon pairing 



e 



nPBF _ « « 1 n28 



6.6 X 102^(m;/m„)(A„(T)/MeV)^ u^/^ 



/(A„(T)/T) erg cm-^-^ (44) 
ePPBF ^ 0.8 X 10^^{m;/mp){Ap{T)/MeVy v"'^ 
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X I{\{T)/T) erg cm-^s'^, 



(45) 



where 

/(A,(r)/r) ~ 0.89 7T/Ai(T)exp[-2A,(T)/T] , (46) 

^ ~ 0.5 for ^Sq pairing and ^ — 1 for pairing. 

A significant contribution of the proton channel is due to the NN corre- 
lation effects, which are taken into account in [H^ . 

The phonon contribution to the emissivity of the 3P2 superfluid phase is 
negligible. The main emissivity regulators are the MMU Eq. (|43|) and the 
neutron (nPBF) and proton (pPBF) pair breaking and formation processes 
Eq. (03) 5 see above for a rough estimation. Finally, we include the effect of 
a pion condensate (PU process), where the PU emissivity is about 1-2 orders 
of magnitude smaller than the DU one. 

All emissivities are corrected by correlation effects. We adopt the same 
set of partial emissivities as in the work of j67j. 

In quark matter the process d — > m + e~ + i^e is always allowed unless it 
is suppressed by huge energy gaps due to quark pairing. 

e«°u ^ 9.4 X lO^^a.^yVs^Qj^^ 7^ ^ni-^ s'S (47) 

where at a compression -u = n/no — 2 the strong coupling constant is ~ 1 
and decreases logarithmically at still higher densities. The nuclear saturation 
density is uq = 0.17 fm~^, = rie/n is the electron fraction, and Tg is the 
temperature in units of 10^ K. If, for a somewhat higher density, the electron 
fraction was too small [Ye < Yec — 10~®), then all the QDU processes would 
be completely switched off [HHI and the neutrino emission would be governed 
by two-quark reactions like the quark modified Urea (QMU) and the quark 
bremsstrahlung (QB) processes dq — > uqev and gig2 <li<l2^i^, respectively. 
The emissivities of the QMU and QB processes have been estimated as jSH] 

Due to the pairing, the emissivities of QDU processes are suppressed by a 
factor Cqdu ~ exp(— Ag/T) and the emissivities of QMU and QB processes 
are suppressed by a factor Cqmu ~ exp(— 2A(,/T) for T < Tent,,; — 0.4 Ag 
whereas for T > Tcrit,q these factors are equal to unity. The modification of 
Tci.it,g(Ag) relative to the standard BCS formula is due to the formation of 
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correlations as, e.g., instanton- anti-instanton molecules [TJ. For the tem- 
perature dependence of the gap below Tcrit,, we use the interpolation formula 
A(T) = A(0)y''l — T/Tcrit,q, with A(0) being the gap at zero temperature. 
The contribution of the reaction ee —>■ eevv is very small [ZO! 

el^ = 2.8 X 10^2 y Vs^i/sj^s ^^-3 g-i^ ^49^ 

but can become important, when quark processes are blocked out for large 
values of A^/T in superconducting quark matter. 



Thermal conductivity 

The heat conductivity of the matter is the sum of the partial contributions 

[zaizs! 

= E'«i' - = E— ' (50) 

where i,j denote the components (particle species). 
The contribution of neutrons and protons is 

= 8.3 X 10^2 erg s-W^K-S (51) 

\m„y bkn -1-9 

Skn = 0.38^-^/2^3.7^2/5 , (52) 
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8.9xl0iM^) erg s- W^K-S (53) 

\m*J ZpSkp 

Skp = 1.83 z-^ + 1.43 4(0.4 + zl)-^ . (54) 

Here we have introduced the appropriate suppression factors Q which 
act in the presence of gaps for superfluid hadronic matter (see Fig. E)) and 
we have used the abbreviation Zi = (nj/(4no))^/^. The heat conductivity of 
electrons is given by Eq. ((HZj). The total contribution related to electrons is 
then 

l/^e = l//tee + l//tep- (55) 

Similar expressions we have for neutrons and protons using Eqs. fl51|) - (jM|) . 

The total thermal conductivity is the straight sum of the partial contri- 
butions Kfot = He + i^n + ■■■ ■ Other contributions to this sum are smaller 
than those presented explicitly (Kg and k„). 
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For quark matter k is the sum of the partial conductivities of the electron, 
quark and gluon components [721 UM 

K = Ke + Kq + Kg, (56) 

where Kg — i^ee is determined by electron-electron scattering processes since 
in superconducting quark matter the partial contribution l/K^q (as well as 
1 / Kgq ) is additionally suppressed by a Cqdu factor, as for the scattering on 
impurities in metallic superconductors. For /tee we have 

Kee = 5.5 X lO^^n Ye Tg'^ crg s"^cm"^K-\ (57) 

and 

Kg - ^ 1.1 X 10^^^^ u (QDuTg-'eig s-^cm-iR-i, (58) 

where we take into account the suppression factor. We estimate the contri- 
bution of massless gluons as 

Kg ~ Kgg ~ 6.0 X 10^^ Big s~^cm"^K"\ (59) 
Heat capacity 

The heat capacity contains nucleon, electron, photon, phonon, and other 
contributions. The main in-medium modification of the nucleon heat capacity 
is due to the density dependence of the effective nucleon mass. We use 
the same expressions as jHIj. The main regulators are the nucleon and the 
electron contributions. For the nucleons {i = n,p), the specific heat is jH] 

c, ~ 1.6 X 102°(m*/m,) (n,/no)'/'C« erg cm'^^K"^ , (60) 

and for the electrons it is 

Ce ~ 6 X 10^^ {ue/nof^^ Tg erg cm"^ K"^ . (61) 

Near the phase transition point the heat capacity acquires a fluctuation con- 
tribution. For the first order pion condensation phase transition this ad- 
ditional contribution contains no singularity, unlike the second order phase 
transition, see [Hni- Finally, the nucleon contribution to the heat capacity 
may increase up to several times in the vicinity of the pion condensation 
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point. The effect of tliis correction on global cooling properties is rather 
unimportant. 

The symmetry of the 3P2 superfluid phase allows a Goldstone boson 
(phonon) contribution of 

cg ^ 6 ■ 10^^T| ergcm-3 K-\ (62) 

for T < Tcn(3P2), n > ncn{<iP2)- We include this contribution in our study 
too, although its effect on the cooling is rather minor. 
For the quark specific heat we use the expression jH^l 

Cq ~ lO^^M^/^Cs Tg erg cm^^ K-\ (63) 

where Cs - 3.1 {T^nt,g/Tf^^ exp{-Ag/T). Besides, one should add the gluon- 
photon contribution 

Cg_y = 3.0 X 10^^ Ng_^ T| erg cm"=^ K-\ (64) 

where Ng_^ is the number of available massless gluon-photon states (which 
are present even in the color superconducting phase), as well as the electron 
contribution, 

4 Results of simulations for Cooling evolu- 
tion 

We compute the neutron star thermal evolution adopting our fully general 
relativistic evolutionary code. The code originally constructed for the de- 
scription of hybrid stars by [T7j has been developed and updated according 
to the modern knowledge of inputs in |18j . 

The density n ~ 0.5 -r- 0.7 no is the boundary of the neutron star interior 
and the inner crust. The latter is constructed of a pasta phase discussed by 
[73] , see also recent works of [7^1 EZI • 

Further on we need a relation between the crust and the surface tem- 
perature for neutron star. A sharp change of the temperature occurs in the 
envelope. This Tg — Tin relation has been calculated in several works, see 
[751 EH] , depending on the assumed value of the magnetic field at the surface 
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and some uncertainties in our knowledge of the structure of the envelope. For 
applications we use three different approximations: the simplified "Tsuruta 
law" T^^^"^ = (10 Tin)^''^ used in many old cooling calculations, models used 
by Ref . |E3 ? and our fit formula interpolating two extreme curves describing 
the borders of region of available Tg — Tin relation, taken from [7H] (see Fig. 
4 of Ref. [Hni)- 

4.1 Cooling Evolution of Hadronic Stars 

Here we will shortly summarize the results on hadronic cooling. 

In framework of a "minimal cooling" scenario, the pair breaking and 
formation (PBF) processes may allow to cover an "intermediate cooling" 
group of data (even if one artificially suppressed medium effects) These 
processes are very efficient for large pairing gaps and temperatures being not 
much less than the value of the gap. 

Gaps, which we have adopted in the framework of the "nuclear medium 
cooling" scenario, see jHO], are presented in Fig. El Thick dashed lines show 
proton gaps which were used in the work of jHj performed in the framework 
of the "standard plus exotics" scenario. We will call the choice of the "3nt" 
model from [Hj the model I. Thin lines show 15*0 proton and 3P2 neutron gaps 
from [HH, for the model AV18 by |H2] (we call it the model II). Recently |^ 
has argued for a strong suppression of the 3P2 neutron gaps, down to values 
~ 10 keV, as the consequence of the medium-induced spin-orbit interaction. 

These findings motivated ^D] to suppress values of 3P2 gaps shown in 
Fig. El by an extra factor f{3P2,n) = 0.1. Further possible suppression of 
the 3P2 gap is almost not reflected on the behavior of the cooling curves. 

Contrary to expectations of jHS| a more recent work of jHll argued that the 
3P2 neutron pairing gap should be dramatically enhanced, as the consequence 
of the strong softening of the pion propagator. According to their estimate, 
the 3P2 neutron pairing gap is as large as 1 10 MeV in a broad region of 
densities, see Fig. 1 of their work. Thus results of calculations of |HS1 and 
[Mj , which both had the same aim to include medium effects in the evaluation 
of the 3P2 neutron gaps, are in a deep discrepancy with each other. 

• Including superfluid gaps we see, in agreement with recent microscopic 
findings of jHS], that the 3P2 neutron gap should be as small as 10 keV 
or less. So the "nuclear medium cooling" scenario of [101 supports 
results of and fails to appropriately fit the neutron star cooling 
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Figure 7: Cooling curves for the gaps from Fig. IHl for model II (the same 
as Fig. 21 of IHOl)- The original 3P2 neutron pairing gap is additionally 
suppressed by a factor /(3P2, n) = 0.1. The Tg — Tin relation is given by "our 
fit" curve of Fig. 4 in jSU]. For more details see [5U] . 



data when a strong enhancement of the 3P2 neutron gaps is assumed 
as shown by |H3] . 

• Medium effects associated with the pion softening are called for by the 
data. As the result of the pion softening the pion condensation may 
occur for n > n^^ {n > Suq in our model). Its appearance at such 
rather high densities does not contradict to the cooling data (see Fig. 
Ej), but also the data are well described using the pion softening but 
without assumption on the pion condensation. The similar argumenta- 
tion holds also for DU threshold density. That puts restrictions on the 
density dependence of the symmetry energy. Both statements might 
be important in the discussion of the heavy ion collision experiments. 

• We demonstrated a regular mass dependence for neutron stars with 
masses M > 1 Mq, where less massive neutron stars cool down slower 
and more massive neutron stars cool faster. This feature is more general 
and concerning also to hybrid stars j86j . 
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Figure 8: Model I. Cooling curves for hybrid star configurations with Gaus- 
sian quark matter core in the 2SC phase. The labels correspond to the 
gravitational masses of the configurations in units of the solar mass. 



4.2 Cooling Evolution of Hybrid Stars with 2SC Quark 
Matter Core 

For the calculation of the cooling of the quark core in the hybrid star we 
use the model (TH]. We incorporate the most efficient processes: the quark 
direct Urea (QDU) processes on unpaired quarks, the quark modified Urea 
(QMU), the quark bremsstrahlung (QB), the electron bremsstrahlung (EB), 
and the massive gluon-photon decay (see ^3]). Following \E7] we include the 
emissivity of the quark pair formation and breaking (QPFB) processes, too. 
The specific heat incorporates the quark contribution, the electron contri- 
bution and the massless and massive gluon-photon contributions. The heat 
conductivity contains quark, electron and gluon terms. 

The calculations are based on the hadronic cooling scenario presented in 
Fig. [7| and we add the contribution of the quark core. For the Gaussian 
form-factor the quark core occurs already for M > 1.214 Mq according to 
the model |17j, see Fig. Most of the relevant neutron star configurations 
(see Fig. [7j) are then affected by the presence of the quark core. 

First we check the possibility of the 2SC+ normal quark phases, see Fig. 

El 
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The variation of the gaps for the strong pairing of quarks within the 2SC 
phase and the gluon-photon mass in the interval 'mg_^ ~ 20 -i- 200 MeV only 
slightly affects the results. The main cooling process is the QDU process on 
normal quarks. We see that the presence of normal quarks entails too fast 
cooling. The data could be explained only, if all the masses lie in a very 
narrow interval (1.21 < M/Mq < 1.22 in our case). In case of the other two 
crust models the resulting picture is similar. 

The existence of only a very narrow mass interval, in which the data can 
be fitted seems to us unrealistic by itself. Moreover the observations show 
existence of neutron stars in binary systems with very different masses, e.g., 
Mbi913+i6 ^ 1.4408 ± 0.0003 Mq and Mjo737_3039b ^ 1-250 ± 0.005 Mq, cf. 
jHHl- Thus the data can't be satisfactorily explained. 

Then we assume a possibility of an yet unknown X-pairing channel with a 
small Ax gap and first check the case Ax to be constant. For Ax — 1 MeV 
the cooling is too slow JH]- This is true for all three crust models. Thus the 
gaps for formerly unpaired quarks should be still smaller in order to obtain 
a satisfactory description of the cooling data. 

For the Ax = 30 keV the cooling data can be fitted but have a very 
fragile dependence on the gravitational mass of the configuration. Namely, 
we see that all data points, except the Vela, CTA 1 and Geminga, correspond 
to hybrid stars with masses in the narrow interval M = 1.21 -^ 1.22 Mq 

Therefore we would like to explore whether a density- dependent X-gap 
could allow a description of the cooling data within a larger interval of com- 
pact star masses. 

We employ an ansatz for the X-gap as a decreasing function of the chem- 
ical potential 

Axifi) = Ac exp[-a(/i - Aic)//ic] , (66) 

where the parameters are chosen such that at the critical quark chemical 
potential /ic = 330 MeV the onset of the deconfinement phase transition 
for the X-gap has its maximal value of Ac = 1.0 MeV and at the highest 
attainable chemical potential /Xmax = 507 MeV, i.e. in the center of the 
maximum mass hybrid star configuration, it falls to a value of the order of 
10 keV. We choose the value a = 10 for which Ax(yUmax) = 4.6 keV. In Fig. 
El we show the resulting cooling curves for the gap model II with gap ansatz 
eq. ()66|). which we consider as the most realistic one. 

We observe that the mass interval for compact stars which obeys the 
cooling data is constrained between M = 1.32 Mq for slow coolers and 
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M = 1.75 Mq for fast coolers such as Vela. This results we obtain are 
based on a purely hadronic model with different choices of the parameters 
|HU] . Note that according to a recently suggested independent test of cooling 
models (HH] by comparing results of a corresponding population synthesis 
model with the Log N - Log S distribution of nearby isolated X-ray sources 
the cooling model I did not pass the test. Thereby it would be interesting to 
see, whether our quark model within the gap ansatz II could pass the Log N 
- Log S test. 




1 2 3 4 5 6 

log(t[yrl) 



Figure 9: Cooling curves for hybrid star configurations with Gaussian quark 
matter core in the 2SC phase with a density dependent pairing gap according 
to Eq. for model II. 



5 Conclusions 

In this lecture devoted to the modern problems of the cooling evolution of 
neutron stars we have discussed different successful scenarios aiming to ex- 
plain the known temperature - age data from observations of compact ob- 
jects. Using up today known theoretical and experimental constraints on the 
structure and cooling regulators of compact star we end up with alternative 
explanations. The neutron star can be either pure hadronic or hybrid one 
with a quark core in superconducting state of matter. 
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Particularly, for the hybrid stars, which are more intrigued alternative of 
superdense compact objects we conclude that: 

• Within a nonlocal, chiral quark model the critical densities for a phase 
transition to color superconducting quark matter can be low enough 
for these phases to occur in compact star configurations with masses 
below 1.3 M0. 

• For the choice of the Gaussian form-factor the 2SC quark matter phase 
arises at M ~ 1.21 Mq. 

• Without a residual pairing the 2SC quark matter phase could describe 

the cooling data only if compact stars had masses in a very narrow 
band around the critical mass for which the quark core can occur. 

• Under assumption that formally unpaired quarks can be paired with 
small gaps Ax < 1 MeV (2SC+X pairing), which values we varied in 
wide limits, only for density dependent gaps the cooling data can be 
appropriately fitted. 

So the present day cooling data could be still explained by hybrid stars, 
assuming a complex pairing pattern, where quarks are partly strongly paired 
within the 2SC channel, and partly weakly paired with gaps Ax < 1 MeV, 
which are rapidly decrease with an increase of the density. 

It remains to be investigated which microscopic pairing pattern could 
fulfill the constraints obtained in this work. Another indirect check of the 
model could be the Log N - Log S test. 
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